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Abstract 

We present results of numerical simulations of the flux (irradiance), F, and 
the degree of polarization (i.e. the ratio of polarized to total flux), P, of light that 
is reflected by Earth-like extrasolar planets orbiting solar-type stars, as functions 
of the wavelength (from 0.3 to 1.0 yum, with 0.001 ^m spectral resolution) and 
as functions of the planetary phase angle. We use different surface coverages for 
our model planets, including vegetation and a Fresnel reflecting ocean, and clear 
and cloudy atmospheres. Our adding-doubling radiative transfer algorithm, which 
fully includes multiple scattering and polarization, handles horizontally homoge- 
neous planets only; we simulate fluxes and polarization of horizontally inhomoge- 
neous planets by weighting results for homogeneous planets. Like the flux, F, the 
degree of polarization, P, of the reflected starlight is shown to depend strongly on 
the phase angle, on the composition and structure of the planetary atmosphere, on 
the reflective properties of the underlying surface, and on the wavelength, in par- 
ticular in wavelength regions with gaseous absorption bands. The sensitivity of P 
to a planet's physical properties appears to be different than that of F. Combining 
flux with polarization observations thus makes for a strong tool for characterizing 
extrasolar planets. The calculated total and polarized fluxes will be made available 
through the CDS. 
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1 Introduction 

Polarimetry has been recognized as a powerful technique for enhancing the con- 
trast between a star and an exoplanet, and hence for the direct detection of exoplan- 
ets, because, integrated over the stellar disk, the direct light of a solar type star can 
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be considered to be unpolarized (see lKemp et al.ll987h . while the starlight that has 
been reflected by a planet will generally be polarized because it has been scattered 
within the planetary atmosphere and/or reflected by the surface (if there is any). 
The degree of polarization of the reflected starlight (i.e. the ratio of the polarized 
to the total flux) is expected to be especially large around a planet's quadrature 
(i.e. when the planet is seen at a phase angle of 90°), where the angular separation 
between a star and its exoplanet is largest, and which is thus an excellent phase 
angle for the direct detection of light from an exoplanet. 

Besides from detecting exoplanets, polarimetry can also be used to character- 
ize exoplanets, because the planet's degree of polarization as a function of wave- 
length and/or planetary phase angle is sensitive to the structure and composition of 
the planetary atmosphere and underlying surface. This application of polarimetry 
is well-known from remote-sen sing of solar system planets, in particular Venus 
(see Hansen & Hovenier 1974a 'b', for a classic example), but also the outer planets 
(see Shkuratov et al. (2005) for recent, H ub ble Space Telescope polarization ob- 
servations of Mars, and ljoos et al. I ( l2005h and lSchmid et"aLl ( l2006h for Earth-based 
polarimetry of Uranus and Neptune). Note that Venus is much more favourable 
to observe with Earth-based polarimetry than the outer solar system planets, be- 
cause as an inner planet, Venus can be observed from small to large phase angles 
(including quadrature), whereas the outer planets are always seen at small phase 
angles, where the observed light is mostly backscatter ed light and degrees of po- 
larization are thus usually small (see lStam et alj ( |2004|) for examples of the phase 
angle dependence of the degree of polarization of starlight reflected by gaseous 
exoplanets). 

The strengths of polarimetry for exoplanet detection and characterization have 

been r e cognized and described b efore, for example by Seager et al. (2000); Saar & Seagetl 
( l2003l) : lHough & Lucasl ( l2003l) : IStaml ( l2003h : IStam et al J ( 120031 12004, 2005), who 
presented numeric ally calculated fluxe s and degrees of polari zatio n of gaseous ex- 
oplane ts. Note that lSeager et alJ ( l2000l)JSaar & Seaged J2003l) , and lHough & Lucaj 
( l2003h concentrate on polarization signals of exoplanets that are spatially ^nresolvable 
from their star, in other words, the polarized flux of the pl anet is added to a huge 
background of unpolarized stellar flux, while StamI ( l2003h and lStam et al.l ( l2003l 
l2004 l 2005h aim at spatially resolvable planets, which are observed with a sig- 
nificantly smaller unpolarized, stellar background signal. Pola rization signals of 
spatial ly unresolved non-spherical planets were presented by ISengupta & Maitil 
( l2006l) . Note that their calculations include only single sca ttered light, and n ot 
all orders of scattering (like those o f ISeager et all ( l2000t) and lStam et al] ( l2004t) ). 
which, except for planetary atmospheres with a very thin scattering or a very thick 
absorption optical thickness, significantly influences the predicted degree of po- 
larization, because multiple scattered light usually has a (much) lower degree of 
polarization than singly scattered light. 

Examples of ground-based telescope instruments that use polarimetry for exo- 
planet research are PlanetPol, whic h aims at detecting spatially unresolved gaseous 
exoplanets (see lHough et alj2006al lbl and references t herein) and SPHER E (Spectro- 
Polarimetric High-contrast Exoplanet Research) (see ISeuzit et al.ll2006l and ref- 
erences therein) which aims at detecting spatially resolved gaseous exoplanets. 
SPHERE is being designed and build for ESO's Very Large Telescope (first light 
is expected in 2010) and ha s a polarimeter based on the ZIMPOL ( Zurich Imaging 
Polarimeter) technique (see lSchmid et al ] |2005l:lGisleretal.ll2004 and references 
therein). Polarimetry is also a technique used in SEE-COAST (the Super Earths 
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Explorer - Coronographic Off- Axis Telescope), a space-based telescope for the de- 
tection and the characteriz ation of gaseous exop lanets and large rocky exoplanets, 
so called 'Super-Earths', jSchneider et alJuOOm . that has been proposed to ESA 
in response to its 2007 Cosmic Vision call. 

The design and development of instruments for the direct detection of (polar- 
ized) light of exoplanets requires sample signals, i.e. total and polarized fluxes as 
functions of the wavelength and a s functions of the planetary phase angle. Pre- 
viously (see e.g. IStam et al]|2004) . we presented numerically calculated flux and 
polarization spectra of light reflected by giant, gaseous exoplanets, integrated over 
the illuminated and visible part of the planetary disk, for various phase angles. 
In this paper, we present similar spectra but now for light reflected by Earth-like 
exoplanets. Our radiative transfer calculations fully include single and multiple 
scattering and polarization. The model atmospheres contain either only gaseous 
molecules or a combination of gas and clouds. The clouds are modeled as hori- 
zontally homogeneous layers of scattering (liquid water) particles, which allows 
surface features to show up in the reflected light even if the planet is fully covered. 
We show results for surfaces with wavelength independent albedos ranging from 
0.0 to 1.0, as well as for surface albedos representative for vegetation, and ocean. 
The ocean surface includes Fresnel reflection. 

Our disk integration method is based on the expansion of the radiation field 
of the planet into generalized spherical functions jStam et al]|2006h . and pertains 
to horizontally homogeneous planets only (the planetary atmospheres can be ver- 
tically inhomogeneous). The main advantage of this method, compared to more 
conventional integration of calculated fluxes and polarization over a planetary disk 
is that the flux and polarization of a planet can be rapidly obtained for an arbitrary 
number of planetary phase angles, without the need of new radiative transfer calcu- 
lations for every new phase angle. This is indeed an important advantage, because 
polarization calculations are generally very computing time consuming compared 
to mere flux calculations. The disadvantage of our method is obviously its inability 
to handle horizontally inhomogeneous planets. In this paper, we will approximate 
the light reflected by horizontally inhomogeneous planets by using weighted sums 
of light reflected by horizontally homogeneous planets. With such quasi horizon- 
tally inhomogeneous planets, we can still get a good impression of the influence of 
horizontal inhomogeneities on the reflected signals. When in the future direct ob- 
servations of Earth-like exoplanets become available, the more conventional disk 
integration method can straightforwardly be applied. 

Our numerical simulations cover the wavelength region from 0.3 to 1.0 yum, 
thus from the UV to the near-infrared. The spectral resolution of our simulations 
is 0.001 yum, which is high enough for spectral features due to absorption of at- 
mospheric gases to be clearly visible in the flux and polarization spectra. Such 
high spectral resolution observations of Earth-like exoplanets will not be possible 
for years to come; our spectra, however, show the potential information content 
of high spectral resolution spectra, and they can be convolved with instrument re- 
sponse functions to simulate observations by instruments with a lower spectral res- 
olution. To allow the use of our flux and polarization spectra for such applications, 
they will be made available at the CDS. 

Flux spe ctra of light reflecte d b y Earth-like exoplanets have been presented be- 
fore ( e.g. bv lTinetti et alj2006bllall3lMontaiies-Rodrfguez et al.l2006l : lTumbull et"aLl 
l2006h . New in this paper are flux spectra with the corresponding polarization spec- 
tra. Numerically calculated degrees and directions of polarization of exoplanets are 
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not only useful for the design and building of polarimeters for exoplanet research, 
as described above, but also for the design, building, and use of instruments that 
aim at measuring only fluxes of exoplanets. Namely, unless carefully corrected 
for, the optical components of such instruments will be sensitive to the state of po- 
larization of the observed light. Consequently, the measured fluxes will depend on 
the state of polarization of the observed light. Provided an instrument's polariza- 
tion sensitivity is known, our simulations can help to estimate the error that arises 
in the measured fluxes. Note that in order to actually correct measured fluxes for 
instrumental polarization sensitivity, knowing the polarization sensitivity of an in- 
strument does not suffice; the state of polarization of the incoming light should be 
measured along with the flux (see e.g. lStam et alj|2000bl . for a discussion on flux 
err ors in remote- sensing due to instrumental polarization sensitivity). 

IStam & HovenieJ | |2005|) discuss another reason to include polarization into 
numerical simulations of light reflected by exoplanets: neglecting polarization in- 
duces errors in numerically calculated fluxes (thus also in e.g. the planet's albedo). 
The reason for these errors is that light can only be fully described by a 4-vector 
(see Sect. 12. U , and a scattering process is only fully described by a 4 x 4 matrix. 
Consequently, the flux resulting from the scattering of unpolarized light differs 
usually from the flux resulting from the scattering of polarized light. Because 
the unpolarized starlight that is incident on a planet is usually polarized upon its 
first scattering, second and higher ord ers of scattering induce errors in the fluxe s 
when polarization is neglected (see also lLacis et al Jj998';' Mi sh chenko et al .Il994l) . 
For gaseous exoplanets, with their optically thick atmosphere s, the flux erro rs due 
to neglecting polarization can reach almost 10 % fstam & HovenieJIlOOSi) . For 
Earth-like exoplanets, with optically thinner atmospheres, we show in this paper 
that the errors are smaller: typically a few percent at short wavelengths ( 0.4 ^m) 
and they decrease with wavelength (see Sect.[5)- 

This paper has the following structure. In Sect.|2] we describe how we define 
and calculate flux vectors and polarization for extrasolar planets. In Sect. |3] we 
describe the atmospheres and surfaces of our Earth-like model extrasolar planets. 
In Sect|4] we present the numerically calculated fluxes and degrees of polarization 
of starlight that is reflected by our Earth-like model planets for both horizontally 
homogeneous planets and the so-called quasi horizontally inhomogeneous planets, 
i.e. weighted mixtures of light reflected by horizontally homogeneous planets. 
Section[5] finally, contains the summary and discussion of our results. 

2 Describing and calculating reflected starlight 
2.1 Flux vectors and polarization 

The flux (irradiance) and state of polarization of stellar light that is reflected by a 
planet can fully be described by a flux (column) vector F as follows 

F = [F, Q, U, V] . (1) 

Here, F is the total reflected flux divided, Q and U des cribe the linearly polar- 
ized flux, and V the circularly polarized flux (see e.g. [Hansen &Travisl [19741 : 
lHovenieretai]|2004h . The fluxes F, Q, U, and V in Eq. [T] have the dimension 
Wm-2 m '. 
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Fluxes Q and U are defined with respect to a reference plane, for which we, 
unless stated otherwise, chose the so-called planetary scattering plane, i.e. the 
plane through the centers of the star and the planet, that also contains the observer. 
We define Fx" as the flux that is measured through a polarization filter oriented 
perpendicular to the direction of propagation of the light, and with its optical axis 
making an angle of x° with the reference plane. The angle is measured rotating 
from the reference plane to the filter's optical axis in the anti-cl ockwise direction 
when looking in the direct ion of propagation of the light (see [Hansen & Travis! 
I1974I : iHovenier et al.ll200 4'). The so-called (linearly) polarized fluxes, Q and U, 
can then in principle be obtained with the following flux measurements 

Q = Foo-F,o", (2) 

U = F45"-Fi35". (3) 

Expressed in the fluxes of Eqs.|2]and[3] the total flux, F, is simply equal to either 
Fo= + Fjo" or F^y + F1350 . Note that modem polarimetry has much more options 
available than polarization filters, such as various types of modulators (see e.g. 
iGandorfer et"ai]|2004 loisler et al.ll2004 ISchmid et al.ll2005l : lKellej2006l) . 

Flux vectors can be transformed from one reference plane to another, e.g. from 
the planetary scattering plane (which depends on the location of the planet on the 
sky with respect to its star) to the optical plane of a polarimeter, by multiplying 
them with a so-called rotation matrix L that is given by (see Hovenier et al. 2004) 



L(yS)^ 



10 

cos2yS sin2y3 

-sin2yS cos2y3 

1 



(4) 



Angle p is the angle between the two reference planes, measured rotating in the 
anti-clockwise direction from the old to the new plane when looking in the direc- 
tion of propagation of the light (j8 > 0). 

The direction of linear polarization with respect to the reference plane is given 
by angle which can be found using 

i2.n2x = UIQ, (5) 

where the convention is to cho ose y such that < y < n, and such that co s 2x 
and Q have the same sign (see [Hansen & Travis|[T974l : iHovenier et alj|2004h . In 
particular, when;^^ = 90° (x = 0°), Q < (Q > 0), and the direction of polarization 
is perpendicular (parallel) to the reference plane, i.e. perpendicular (parallel) to 
the imaginary line connecting the centers of the star and the planet as seen from 
the observer. 

The degree of polarization of the reflected starlight is defined as 

P^Ml^. (6) 
F ^ ' 

Note that P as defined in Eq. [6] is independent of the choice of reference plane. 
Assuming that the planet is mirror-symmetric with respect to the planetary scatter- 
ing plane, and assuming the incoming starlight is unpolarized, the disk-integrated 
Stokes parameters U and V will equal zero because of symmetry (the incoming 
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starlight is unpolarized). In that case, we can use an alternative definition of the 
degree of polarization (see also Eq.|2j, namely 

P^ = -^ = -^L^, (7) 

t too + tgo- 

For Ps > 0, the light is polarized perpendicular to the reference plane (i.e. x = 90°), 
and for > 0, the light is polarized parallel to the reference plane (i.e. ^ = 0°). 
The subscript i (from "sign") in P^ thus indicates that the direction of polarization 
is included in the definition of the degree of polarization. 



2.2 Computing reflected starlight 

Given a spatially unresolved spherical planet with radius r, the flux vector F (see 
Eq. [T]l of stellar light with wavelength A that has been reflected by the planet and 
that arrives at an observer at a distance d (with d » r) can be written as (see also 
IStam et al.ll2006h 

r2 1 

F(^,a)= - -S(^,ff);rFo(^). (8) 
d- 4 

Here, a is the planetary phase angle, i.e. the angle between the star and the observer 
as seen from the center of the planet. Note that a = 180° - 0, with the total 
scattering angle of the incoming starlight. A sketch of the geometries is given in 

Fig.m 

Furthermore in E q. [S] S is the 4x4 planetary scattering matrix (see below 
and IStam et alj|2006 1 and Fo represents the flux (column) vector describing the 
stellar light that is incident on the planet, with nFo the stellar flux that arrives at 
the planet (in W m"- m ') measured perpendicular to the direction of propagation 
of the stellar light. Integrated over the stellar disk, the stellar light of a solar-type 
star can be assumed to be unpolarized (Kemp et al. 1987), hence in the following 
we will use Fo(A) = Fo(/l)[l,0, 0, 0] = Fo(A)l, with 1 the unit column vector. We 
assume that the starlight is unidirectional when it arrives at the planet. 

The planetary scattering matrix, S, depends on the planetary phase angle, a, 
and on the wavelength, A. The relation between S and a and A depends on the com- 
position and structure of the planetary atmosphere and on the planetary surface. 
This dependence will be further described in Sect. [3] Using the planetary scatter- 
ing plane as the reference plane, and assuming the planet is mirror-symmetric w ith 
respect to this reference plane, matrix S is given by (see lStam et al.ll200^l2004b 



S(A,a) : 



a,(A,a) bi(A,a) 

bi(A,a) a2{A,a) 

ai(A,a) b2(A,a) 

-h2(A,a) a4(A,a) 



(9) 



Matrix element ai is usually called the planetary phase function. Matrix S is nor- 
malized such that the average of aj over all directions equals the planet's (monochro- 
matic) Bond albedo, Ab, which is the fraction of the incident stellar flux that is 
reflected by the planet in all directions, i.e. 

If I r" 

— I ai(A,a) dcij = — I ai(A,a) sin a da = A^(A), (10) 
4?r Ja^ 2 Jo 

where da) is an element of solid angle. The (monochromatic) geometric albedo, 
Aq, of a planet is the ratio of the flux reflected by the planet at a = 0°, to the flux 
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reflected by a Lambertian surface subtending the same solid angle (i.e. nr^ 1(f) on 
the sky. Thus, 

F(/l,0°) £- 1 

^o('i> = 4-hr ^ = 7'''('^'0 )• (11) 

Fo(/i) nr'- 4 

From Eqs. [8]and|9]it is clear that with unpolarized incident stellar light, the 
observable total flux, F, of starlight that is reflected by a planet is given by 

F(A,a)= ^ LiiA,a)nFo(A), (12) 
and the observable polarized flux, Q, by 

Q(A,a) = Ui(A,a)nFo(A). (13) 

With unpolarized incident stellar light and a planet that is mirror-symmetric 
with respect to the planetary scattering plane, it follows from Eas.ll2landll3lthat 
the degree of polarization of the starlight that is reflected by the planet can 
simply be rewritten as (cf. Eq.|7]l 

P.(,,.) = -^. (14) 

The degree of polarization of the reflected light then thus solely depends on the 
planetary scattering matrix elements and bi. Because both P (Eq.[6} and P^ 
(Eqs.[7landll4t are relative measures, they are independent of the radii r and R, the 
distances D and d, which is very convenient when analyzing direct observations of 
extrasolar planets at unknown distances. 

To calculate the flux and degree of polarization of light reflected by a given 
Earth-like model planet (see Sect. [3]l across a given wavelength region and for a 
given planetary phase angle, we have to calculate elements of the planetary s cat- 
tering matrix S (Eq. |9). For this we use the algorithm as descr ibed i n Stam et al.| 
(I2OO6I), which combines an accurate adding-doubling algorithm ( ivan de Hulstll980l ; 
ldeHaanetaLlll987l) to compute the radiative transfer through a locally plane- 
parallel planetary model atmosphere, and a fast, numerical, disk-integration al- 
gorithm, to integrate the reflected flux vectors across the illuminated and visible 
part of the planetary disk. 

Our disk-integration algorithm i lStam et alj20o3) is very efficient, and its com- 
puting time depends only little on the number of planetary phase angles for which 
the disk-integrated flux vectors are calculated. The disadvantage of the current 
version of the algorithm is that it can only handle horizontally homogeneous plan- 
ets (which are mirror-symmetric with respect to the reference plane). Thus, while 
a planetary model atmosphere can be inhomogeneous in the vertical direction, it 
varies neither with latitude nor with longitude. Calculated flux and polarization 
spectra of horizontally homogeneous planets are presented in Sects. [43] and |4.2| 
In this paper, we will approximate horizontally inhomogeneous planets by using 
weighted sums of horizontally homogeneous planets. The flux vector of such a 
quasi horizontally inhomogeneous planet is calculated according to 

IV N 

F(A,a) = Y,fnF„(A,a) with ^Z, = 1, (15) 

n=l 11=1 

with A' the number of horizontally homogeneous planets. Calculated flux and po- 
larization spectra of quasi horizontally inhomogeneous planets are presented in 
Sect. 143] 
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2.3 Atmospheric extinction and instrumental response 

The flux vector F as described in Eq.[8]includes neither extinction in the terrestrial 
atmosphere, nor the response of an instrument. It thus pertains to the flux vector 
as it can be observed in space. Adding atmospheric extinction and/or instrumental 
effects is straightforward by multiplying vector F from Eq. [8]with the matrix de- 
scribing atmospheric extinction, and/or with the matrix describing the instrumental 
response. 

Atmospheric extinction usually affects only the flux of the directly transmitted 
light, not its state of polarization, and can then simply be described by the scalar 
exp(-T) cos"'(0o), with r the (wavelength dependent) extinction optical thickness 
of the atmosphere between the observer and space, and Oq the zenith angle of the 
observed exoplanet. 

Since most instruments not only affect the fl ux of the light that enters the instru- 
ment, but also its state of polarization, (see e.g. lStam et al.l2000bl . regarding polar- 
ization sensitive Earth-observation instruments) an instrumental response matrix 
can be quite complicated. With a polarization sensitive instrument, the flux that 
is measured will not only depend on the flux F (see Eq. I12t of the light that is 
reflected by the planet, but also on its state of polarization. Of course, a polar- 
ization sensitive instrument will usually also change the state of polarization of 
the observed light. Thus, when analyzing flux and/or polarization observations 
of starlight that is reflected by an exoplanet, one has to properly account for the 
polarization sensitivity of one's instrument. 

In this paper, we will ignore both atmospheric extinction and instrumental ef- 
fects (apart from the spectral resolution of 0.001 fim), and thus limit ourselves to 
the total flux and degree of polarization as they can be observed in space. Since 
the calculated total and polarized fluxes will be made available through the CDS, 
the atmospheric extinction and/or instrumental effects that particular observations 
require can be applied. 

3 The model planets 

The atmospheres of our Earth-like model planets are described by stacks of ho- 
mogeneous layers containing gaseous molecules and, optionally, cloud particles. 
Each model atmosphere is bounded below by a flat, homogeneous surface. In the 
next subsections, we will describe the composition, structure, and optical prop- 
erties of our model atmospheres (Sect. 13. IK and the reflection properties of our 
model surfaces (Sect. [T2t . 

3.1 The model atmospheres 

All of our model atmospheres consist of 16 homogeneous layers. For the radia- 
tive transfer calculations, we need to know for each atmospheric layer: its optical 
thickness, b, and the single scattering albedo, a, and scattering matrix, Fj^a, (see 
Hovenier et al., 2004) of the mixture of molecules and cloud particles. 

An atmospheric layer's optical thickness, b, is the sum of its molecular and 
cloud extinction optical thicknesses, b™ and b'^, i.e. 

b(A) = b'^HA) + b\X) = bZM) + b2M) + bl-M) + *ab.W- (16) 
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Here, and b^^^^ are the molecular scattering and absorption optical thicknesses, 
respectively, and bl^.^ and b':^^^ are the cloud scattering and absorption optical thick- 
nesses. 

The molecular scattering optical thickness of each atmospheric layer, is 
calculated as described in Stam et al. (2000a), and depends a.o. on the molecular 
column de nsity (molecules per m ^), the refractive index of dry air under standard 
conditions jPeck & Reede3ll972 l). and the depolarizati on factor of a ir, for which 
we adopt the (wavelength dependent) values provided bv lBatesl ( ll984l) . The molec- 
ular column density depends on the ambient pr essure and temp erature, the vertical 
profile of which is given in Table [T] (McClatchev et alj|l972h (to avoid introduc- 
ing too many variables, we use this mid-latitude summer vertical profile for each 
model atmosphere). In Table[T] we also give for each atmospheric layer Z?™^ as we 
calculated at A = 0.55 /jm. 

The molecular absorption optical thickness of each atmospheric layer, fc'^^, de- 
pends on the molecular column density, the mixing ratios of the absor bing gases, 
and their mo lecular absorption cross-section (in m^ per molecule) (see lStam et al.l 
Il999l , [2000al , for the details). The terrestrial atmosphere contains numerous types 
of absorbing gases. In the wavelength region of our interest, i.e. between 0.3 yum 
and 1.0 yum, the main gaseous absorbers (and the only absorbers we take into ac- 
count here) are ozone (O3), oxygen (O2) and water (H2O). Unless stated other- 
wise, the mixing ratio of O2 is 2.1-10'* ppm (parts per million) throughout each 
model atmosphere. The altitude dependent mixing ratios of the trace gases O3 and 
H2O are given in Table [T] We calculate the molecular absorption cross-sections 
of O2, O3, and H2O using absorption line data from iRothman et al.l l l200j) . Be- 
cause the absorption cross-sections of O2 and H2O are rapidly varying functions 
of the wavelength, we have transforme d them into so-called /^-distributions (see 
iLacis & Oin"aslll99ll :[ Stam et aill2000ah . using a wavelength spacing of 0.001 yum, 
a spectral resolution of 0.001 yum, 20 Gaussian abscissae per wavelength interval 
of 0.001 fim, and a block-shaped instrumental response function. The absorption 
cross-sections of O3 vary only gradually with wavelength (at least between 0.3 yum 
and 1.0 /im); we assume them to be constant across each wavelength interval of 
0.001 fim. Molecular absorption cross-sections in general not only depend on the 
wavelength, but also on the ambient pressure and temperature. For the purpose 
of this paper, i.e. presenting flux and polarization spectra of Earth-like extrasolar 
planets and addressing the occurence of spectral features in them, we use the ab- 
sorption cross-sections calculated for the lowest atmospheric layer (see Table [TJ 
throughout our model atmospheres. 

For each wavelength and each atmospheric layer, the cloud scattering and ab- 
sorption optical thicknesses, bl^^ and fo^^^^, are calculated from the user-defined 
cloud particle column density (in cloud particles per m^), and the extinction cross- 
section and the single scattering albedo of the cloud particles. The only cloud 
particles we will consider in this paper are spherical, homogeneous, watercloud 
droplets. These dropl ets are distributed in size according to the standard size dis- 
tribution described by [Hansen & TravisI jl974) . with an effective radius of 2.0 yum 
and an effective variance of 0.1. The refractive index is chosen to be wavelength 
independent and equal to 1.33 + 0.0001. We calculate the extinction cross-section, 
single scattering albedo, and the scattering matrix, F^^.^, of th e cloud droplets for 
wavelengths between 0.3 a nd 1.0 /jm using Mie-theory (see lvan deHulstll957l : 
Ide Rooii & van der Staplfigsj) . Unless stated otherwise, we assume a cloud with 
an optical thickness, b", of 10 at /I = 0.55 fim, with its bottom at 802 hPa and its 
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top at 628 hPa (according to Table[T] it thus extends vertically from 2 km to 4 km). 
Because of the wavelength dependence of the droplets' extinction cross-section, 
the cloud's optical thickness varies with wavelength. In particular, al A = 0.3 yum, 
6" = 9.6, and at ^ = 1.0 yum, // = 10.7. 

The single scattering albedo of the mixture of gaseous molecules and cloud 
particles in an atmospheric layer is calculated according to 

(A) = fcSaW + ^LaW .,7. 

and the scattering matrix (see lHovenier et alj2004) of the mixture according to 

F re A) - + bl/A)FlS0,A) 

' bZM) + blJX) 

where is the scattering angle (with = 0° indicating forward scattering), and 
F™;, and F^^.^ are the scattering matrices of, respectively, the molecules and the 
cloud particl es. The scattering m atrix F™^ of the gaseous molecules is calculated as 



described bv lStam et al (2003), using the (wavelength dependent) depolarization 
factor of air (' Bates! 1984 ). We do not explicitly account for rotational Raman scat- 
tering, an inelastic molecular scattering process (see e.g . jGrainger & Ring 1 963 : 
lAben et al.l200ll : IStam et alj2002l : lvMi Deelen et alj200 j : ISromovskvf200^ . which 
gives rise to a slight "fiUing-in" of high-spectral resolution features in reflected 
light spectra, such as stellar Fraunhofer lines and gaseous absorption bands. Each 
scattering matrix is normalized such that the average of the phase function, which 
is repre sented by scattering matrix element F", ove r all scattering directions is 
one (seeH ansen & Travislll974lHovenier et aU2004h . 

Figure |2^ shows the phase functions of the gaseous molecules and the cloud 
droplets aX A = 0.55 fim. To illustrate the wavelength dependence of the ele- 
ments of the cloud droplets' scattering matrix, we have also plotted curves for 
A = 0.44 /jm and A = 0.87 yum (these particular wavelengths will be used again 
later on, in Sect.|4l(. For the same wavelengths, Fig.|2j) shows the degree of linear 
polarization, (Eq. |7}, of light that is singly scattered by the molecules and the 
cloud droplets as functions of the single scattering angle, 0, assuming unpolar- 
ized incident light. The reference plane for this singly scattered light is the plane 
through the incoming and the scattered light beams. Note that the phase function 
and degree of polarization of light singly scattered by gaseous molecules depends 
on the wavelength, too, through the wavelength dependence of the depolarization 
factor (see Bates 1984), but only slightly so. 

From Fig. |2l it is clear that both the phase function and the degree of polar- 
ization pertaining to single scattering by molecules vary smoothly with the single 
scattering angle 0. The degree of polarization, f s, of the light that is singly scat- 
tered by the molecules is positive (i.e. the direction of polarization is perpendicular 
to the reference plane) for all values of 0. Furthermore, of this light is highest 
at = 90°. At this scattering angle, the light is not completely (i.e. 100 %) polar- 
ized, but "only" about 95 %, because of the molecular depolarization factor ( iBatesI 
Il984h . Both for the light scattered by the molecules and for the light scattered by 
the cloud droplets, f s vanishes for = 0° and = 180°, because of symmetry. 

The degree of polarization of the light scattered by the cloud droplets (Fig.|2]3) 
changes sign (i.e. the direction of polarization changes with respect to the reference 
plane) a number of times between = 0° and 0=1 80° , and shows strong angular 
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features, in particular in the backward scattering directions (0 > 90°). The peak 
in the polarization occuring at = 148° for A = 0.44 fim, at = 150° for 
A = 0.55 ftm, and at = 155° for A = 0.87 /im, pertains to what is commonly 
known as the primary rainbow, which is due to light that has been reflected inside 
the droplets once. The primary rainbow is seen in the flux phase functions (Fig.[2^), 
too, only less prominent than in P^. The angular features in near © = 120° 
pertain to the secondary rainbow, which is due to light that has been reflected inside 
the the droplets twice. In the cloud droplets' phase functions (Fig.|2^), only a hint 
of the secondary rainbow can be seen and only for A = 0.44 fim. The occurence 
of a rainbow in reflected light is a str ong and well-known indic ator for spherically 
shaped atmospheric pa rticles, see e.g. [Hansen & Travis! (Il974h, and more recently 
iLiou & Takand j2002l) and references therein, and BailevI j2007h . 



3.2 The model surfaces 

To describe the reflection of light by the homogeneous, locally flat surfaces below 
the atmospheres of our model planets, we have to specify the surface reflection 
matrix, Aj. The surface reflection matrix is normalized such that the average of 
matrix element (1,1) of As over all reflection directions equals the surface albedo, 
i.e. the fraction of the incident stellar flux that the surface reflects in all directions. 
We'll denote the surface albedo by A^. 

The Earth's surface is covered by numerous surface types with myriads of 
(wavelength dependent) albedos, many of which vary with e.g. their moistness 
and/or the season. To avoid making our model planets too detailed at this stage, 
we will compose the surfaces of our Earth-like model planets out of only two sur- 
face types: (deep) ocean and (green) vegetation. We assume that the surface that is 
covered by vegetation completely depolarizes all incident light, i.e. except for ele- 
ment element (1,1) of As, all elements of the vegetation's surface reflection matrix 
equal zero. In addition, we assume that the reflection by the surface is isotropic, 
i.e. reflection matrix element (1,1) is independent of the directions of both the in- 
coming and the reflected light; element (1,1) thus simply equals the surface albedo. 
As. We thus describe a surface that is covered with vegetation as a Lambertian re- 
flecting surface. In future stud ies, it will be interesting to include polarizing effects 
of vegetation, as presented bv lWolstencroft et alj ( l2007h . 

In Fig. |3] we have plotted measured, wavelength dependent albedos of three 
types of vegetation: conifers, deciduous forest, and grass Q These albedo spec- 
tra share the following characteristics: (1) a local maximum between 0.5 yum and 
0.6 yum, that is mainly due to the presence of two absorption bands of chlorophyll, 
one near 0.45 fim and one near 0.67 fim, and (2) a high albedo at wavelengths 
longer than about 0.7 yum, that is related to the internal leaf and cell structure. The 
sudden increase of the surface albedo at wavelengths longer than 0.7 yum, is usu- 
ally referred to as the red edge (for an elaborate description of the red edge, see 
ISeager et al.ll2005l) . The slight decrease in the albedo around 0.97 fim is due to 
absorption by water in the leaves. Stronger absorption bands of water occur at 
wavelengths longer than 1.2 yum. Because in this paper we do not study the ef- 
fects of differences in the albedos of different types of vegetation on the light that 
is reflected by a planet, we will only use the wavelength dependent albedo of the 

' These three albedos have been taken from the ASTER Spectral Library through the courtesy of the Jet 
Propulsion Laboratory, California Institute of Technology, Pasadena, California. 
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deciduous forest to represent the reflectivity by vegetation on our model planets. 

Whether vegetation on Earth-like extrasolar planets will have the same spec- 
tral featu res, in particular the red edge , as we find on Earth, is still an open ques- 
tion (see IWolstencroft & Ravenll2002h . Model studies for albedos of vegetation 
on Earth-like plane ts around M stars have been published by Kiang et al. (200^; 
ISegura et alj ilOO^ : iTinetti et alj |2006c). For our purpose, presenting flux and 
polarization spectra and in particular their differences and similarities without fo- 
cussing on the detection of features, an Earth-like vegetation albedo is sufficient. 

Although on Earth, deep oceans do show some color, especially in shallow re- 
gions where algae and other small organisms bloom, for the purpose of this paper 
it is safe to simply assume the oceans are black across the wavelength interval of 
our interest, i.e. from 0.3 to 1.0 /jm. Even with an albedo As equal to zero, how- 
ever, our model oceans do reflect a fraction of the light that is incident on them, 
because we include a specular (i.e. Fresnel) reflecting interface between the atmo- 
sphere and the black ocean. Specular reflection is anisotropic and generally leads 
to polarized reflected light. We use the specular reflection matrix as described by 
iHaferman et"^ ( fl997i) . with a (wavelength independent) index of refraction that 
is equal to 1.34. Our model ocean surface is flat, i.e. there are no waves. The 
influence of oce anic waves wil l be su bject of later studies, using the wave distri- 
bution model by ICox & Muri3 ( 1 1954) (for a recent evaluation of this model, see 
iBreon & HenrQl2006h . which can be included in our adding-doubling radiative 
transfer model (see e.g. 'Chowdharv et al.l l2002l) . We neglect the contribution of 
whitecaps to the ocean albedo, which appears to be a valid a ssumption for average 
wind speeds measured on the Earth's oceans ( lKoeDkdiT984b . 

4 Calculated flux and polarization spectra 

In this section, we will present the numerically calculated total flux and degree of 
polarization of starlight that is reflected by Earth-like model planets as described 
in the previous section (Sect. [3}. The reflected flux, F, is calculated according 
to Eq. [8] Unless stated otherwise, we assume that r = 1, rf = 1, and nFo = 1, 
independent of A. With unpolarized incident light, F thus equals which is the 
planet's geometric albedo Aq in case planetary phase angle a = 0° (see Eq. II U . 
The degree of polarization is calculated according to Eq.[T4l and thus includes the 
direction of polarization. 

Tables containing elements and bi of the planetary scattering matrix S as 
functions of the wavelength, and as functions of the planetary phase angle, for the 
various horizontally homogeneous model planets that are presented in the follow- 
ing sections will be made available through the CDS. From the elements oi and 
bi, and given distance d, planetary radius r, and the (wavelength dependent) inci- 
dent stellar flux (e.g. in W nr^ m"'), the observable total flux, F, polarized flux, 
Q, and degree of polarization, P^, can be calculated using Eqs. [T2l[l3] and 1 141 
respectively. 
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4.1 Clear planets with wavelength independent surface 
albedos 



4.1.1 Wavelength dependence 

Figure |4] shows the wavelength dependence of the total flux, F, and the degree of 
polarization, P^, of starlight that is reflected by six Earth-like model planets with 
similar, clear (i.e. cloudless) atmospheres, and Lambertian reflecting (i.e. isotrop- 
ically reflecting and completely depolarizing) surfaces with wavelength indepen- 
dent albedos, A^, ranging from 0.0 to 1.0. The planetary phase angle, a, is 90°, 
i.e. half of the observable planetary disk is illuminated by the star. As explained in 
Sect.[T] the probability to directly observe an exoplanet at or near this phase angle 
(quadrature) is relatively high (provided there is an observable exoplanet). 

Each curve in Fig.|4]can be thought of to consist of a continuum with superim- 
posed high-spectral resolution features. The continua of the flux and polarization 
curves are determined by the scattering of light by gaseous molecules in the at- 
mosphere and by the surface albedo. The high-spectral resolution features are 
due to the absorption of light by the gases O3, O2, and H2O (see below). Note 
that strength and shape of the absorption bands depend on the spectral resolution 
(0.001 /jm) of the numerical calculations. 

In the total flux curves (Fig. |4^), the contribution of light scattered by atmo- 
spheric molecules is largest around 0.34 /jm; at shorter wavelengths, light is ab- 
sorbed by O3 in the so-called Huggins absorption band, and at longer wavelengths, 
the amount of starlight that is scattered by the atmospheric molecules decreases, 
simply because the atmospheric molecular scattering optical thickness decrease s 
with wavelength, as i"., is roughly proportional to /l""* (see e.g. IStam et al.l2000ah . 
For the planet with the black surface (As = 0.0), where the only light that is re- 
flected by the planet comes from scattering by atmospheric molecules, the flux 
of reflected starlight decreases towards zero with increasing wavelength. For the 
planets with reflecting surfaces, the contribution of light that is reflected by the sur- 
face to the total reflected flux increases with increasing wavelength. Because the 
surface albedos are wavelength independent, the continua of the reflected fluxes 
become wavelength independent, too, at the longest wavelengths. This is not ob- 
vious from Fig.|4^, because of the presence of high-spectral resolution features. 

The high-spectral resolution features in the flux curves of Fig.|4]are all due to 
gaseous absorption bands. As mentioned above, at the shortest wavelengths, light 
is absorbed by O3 . The so-called Chappuis absorption band of O3 gives a shallow 
depression in the flux curves that is visible between about 0.5 yum and 0.7 yum, 
in particular in the curves pertaining to a high surface albedo. The flux curves 
contain four absorption bands of O2, i.e. the y-band around 0.63 yum, the B-band 
around 0.69 yum, the conspicuous A-band around 0.76 yum, and a weak band around 
0.86 /jm. These absorption bands, except the A-band, are difficult to identify from 
Fig.E^, because they are located either next to or within one of the many absorption 
bands of H2O (which are all the bands not mentioned previously). 

The polarization curves (Fig.|4j5) are, like the flux curves, shaped by light scat- 
tering and absorption by atmospheric molecules, and by the surface reflection. 
The contribution of the scattering by atmospheric molecules is most obvious for 
the planet with the black surface (As = 0.0), where there is no contribution of the 
surface to the reflected light. For this model planet and phase angle, has a local 
minimum around 0.32 yum. At shorter wavelengths, P^ is relatively high because 
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there the absorption of light in the Huggins band of O3 decreases the amount of 
multiple scattered light, which usually has a lower degree of polarization than the 
singly scattered light. In general, with increasing atmospheric absorption optical 
thickness, will tend towards the degree of polarization of light singly scattered 
by the atmospheric constituents (for these model planets: only gaseous molecules), 
which depends strongly on the single scattering angle and thus on the planetary 
phase angle a. From Fig.[2j3, it can be seen that at a scattering angle of 90°, Ps of 
light singly scattered by gaseous molecules is about 0.95. This explains the high 
values of P^ at the shortest wavelengths in Fig.|4]3. With increasing wavelength, the 
amount of multiple scattered light decreases, simply because of the decrease of the 
atmospheric molecular scattering optical thickness. Consequently, P^ of the planet 
with the black surface increases with wavelength, to approach its single scattering 
value at the smallest scattering optical thicknesses. 

With a reflecting surface below the atmosphere, P^ also tends to its single scat- 
tering value at the shortest wavelengths, because with increasing atmospheric ab- 
sorption optical thickness, the contribution of photons that have been reflected by 
the depolarizing surface to the total number of reflected photons decreases (both 
because with absorption in the atmosphere, less photons reach the surface and less 
photons that have been re flected by the surface reach the top of the atmosphere) 
(see e.g. lStamet J]|l999h . In case the planetary surface is reflecting, P^ of the 
planet will start to decrease with wavelength, as soon as the contribution of pho- 
tons that have been reflected by the depolarizing surface to the total number of 
reflected photons becomes significant. As can been seen in Fig. |4]3, the wave- 
length at which the decrease of P^ starts depends on the surface albedo: the higher 
the albedo, the shorter this wavelength. It is also obvious that with increasing 
wavelength, the sensitivity of P^ to decreases. This sensitivity clearly depends 
on the atmospheric molecular scattering optical thickness. 

Like with the flux curves, the high-spectral resolution features in the polariza- 
tion curves of Fig. |4j3 are all due to gaseous absorption. The explanation for the 
increased degree of polarization inside the O2 and H2O absorption bands is the 
same as that given above for the Huggins absorption band of O3: with increas- 
ing atmospheric absorption optical thickness, the contribution of multiple scat- 
tered light to the reflected light decreases, and hence P^ increases towards the 
degree of polarization of light singly scattered by the atmospheric constituents, 
i.e. gaseous molecules. In case atmospheres contain aerosol and/or cloud parti- 
cles, Ps both inside and outside the absorption bands will depend o n the single 



scattering properties of those aerosol and/or cloud particles, too (see Stam etal 
1 19991. for a detailed des cription of P^ across gaseous absorption lines). ,Stam et al. 



j2004h and [staii]b003h show calculated polarization spectra of Jupiter-like extra- 
solar planets with gaseous absorption bands due to methane. An increase of the 
degree of polarization across gaseous asorption bands has been measured in so- 
called zenith sky observations on Earth jStammes et al.l 19941 : IPreusker et aljl995l : 
lAbenetal.ll 19971.1 1999I), an d, recently, in observations of Jupiter , Uranus and Nep- 
tune jjoos & Schmi j|2007l : ISchmid et alj|2006l : ljoos et al.ll2005h . with methane as 
the absorbing gas. 

It is interesting to note that the polarization spectrum of an extrasolar planet 
will generally be msensitive to absorption that takes place between the planet and 
the observer because it is a relative measure (see Eqs.[6]and|7}. Thus, if the tele- 
scope were located on the Earth's surface, the polarization features as shown in 
Fig. |4j) would be unaffected by absorption within the Earth's atmosphere; po- 
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larimetry would in principle allow the detection of e.g. O2 in an extrasolar plane- 
tary atmosphere despite the O2 in the Earth's atmosphere (the number of photons 
received by the telescope, i.e. the flux, would of course be strongly affected by 
absorption in the Earth's atmosphere). 

4.1.2 Phase angle dependence 

Figure |5] shows the phase angle dependence of the total flux, F, and the degree of 
polarization, P^, of the starlight that is reflected by three of the six Earth-like plan- 
ets appearing in the previous section, namely, the planets with As = 0.0, 0.4, and 
1.0. The phase angle dependence has been plotted for two wavelengths: 0.44 yum 
and 0.87 fim. Remember that the flux at phase angle or = 0° is just the planet's 
geometric albedo Aq (Eq. [TTJ- For = 0.0, Aq = 0.14 at /I = 0.44 yum, and 
Ag = 0.011 atA = 0.87 pm. For A^ = 0.4 and 1.0, we find, respectively, Aq = 0.34 
(0.44 yum) and 0.27 (0.87 pm), and 0.72 (0.44 yum) and 0.67 (0.87 yum) (see Fig.[5li). 
The planets' geometric albedos at /I = 0.87 ytim are close to | A,, i.e. the geometric 
albe do of a planet wi th a Lambertian reflecting surface but without an atmosphere 
(see lStam et"ai]|2006h . because at this wavelength, the scattering optical thickness 
of the model atmosphere is only 0.015. At 0.44 yum, this optical thickness is 0.24, 
and the light scattered within the model atmosphere does contribute significantly 
to the planet's geometric albedo, especially when As is small. 

The strong phase angle dependence of F (Fig. |5ji) is, for a given value of As, 
largely due to the variation of the illuminated and visible fraction of the planetary 
disk with the phase angle. Other variations are related to the reflection properties 
of the surface and the scattering properties of the overlying atmosphere. These 
variations can be seen more clearly in Fig. (6) where we show the curves of Fig. [5^ 
normalized at a = 0°. The curves for A = 0.87 yum and As = 0.4 and As = 1.0 co- 
incide with the theoretical (normalized) curve ex pected for Lambertian reflecting 
spheres (see lvan de Hulstiri98(ll : IStam et al.ll200d) . because, as explained above, at 
this wavelength, the contribution of light scattered within the model atmosphere is 
almost negligible. 

In Fig.[5}), we show the degree of polarization, P^, as a function of the phase 
angle. Like the reflected flux, P^ depends strongly on the phase angle. Note that 
for a = 0° and a = 180°, Ps equals zero because of symmetry (the incoming 
starlight is unpolarized). For the planet with the black surface, P^ appears to be 
fairly symmetric around a = 90°, mainly because the degree of polarization of 
light singly scattered by gaseous molecules is symmetric around = 90° (see 
Fig- lib)- This symmetry is particularly apparent for A = 0.87 yum, where there is 
much less multiple scattering than for A = 0.44 yum. Across most of the phase 
angle range, f s of the planet with the black surface is positive, indicating that the 
reflected light is polarized perpendicular to the scattering plane (i.e. perpendicular 
to the imaginary line connecting the planet and the star as seen by the observer). 
Only for the largest scattering angles, P^ < (which is difficult to see in Fig.|5j) 
for A = 0.87 fjm). At these angles, the reflected light is thus polarized parallel 
to the scattering plane (i.e. parallel to the imaginary line connecting the planet 
and the star). The negative polarization is mainly due to second order scattered 
light: although this light comprises only a small fraction of the reflected light at 
these large phase angles, it does leave its traces in the polarization signature of the 
planet, because the degree of polarization of the first order scattered light (which is 
the main contributor to the reflected light) is close to zero. For A = 0.44 yum, Ps < 
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when a > 164°, while for /I = 0.87 yum, Ps < only when a > 174°, because with 
increasing wavelength, the amount of second order scattered light decreases, and 
with that the phase angle at which the second order scattered light changes the sign 
of Ps increases. 

For the planets with the reflecting surfaces, the maximum degree of polar- 
ization occurs at phase angles larger than 90° (see Fig.jSjj). In particular, with 
increasing wavelength and/or increasing surface albedo, the maximum degree of 
polarization shifts to larger phase angles, because with increasing a, the fraction 
of reflected light that has touched the depolarizing surface at least once, decreases. 
The contribution of light that has been polarized within the planetary atmosphere 
to the reflected light thus increases with increasing a. This also explains why at 
the largest values of a, Ps in Fig.[5j3, is independent of A^. 

Figure [7] shows F and Ps as functions of the planet's orbital position angle for 
orbital inclination angles, /, ranging from 0° (the orbit is seen face-on) to 90° (the 
orbit is seen edge-on). Given an inclination angle an exoplanet can in princi- 
ple be observed at phase angles, a, ranging from 90° - ; to 90° + /. If the orbital 
position angle equals 0° or 360°, a ranges from 90° (; = 0°) to 0° (/ = 90°). If 
the orbital position angle equals 90° or 270°, the planetary phase angle, a, equals 
90° (independent of ;)• If the orbital position angle equals 180°, a ranges from 
90° (/ = 0°) to 180° (/ = 90°). The curves in Fig. |7] clearly show that with in- 
creasing orbital inclination angle, the variation of F and Ps along the planetary 
orbit increases. Interestingly, the orbital position angles where a planet is easi- 
est to observe directly because it is furthest from its star (in angular distance) are 
those where Ps is largest (namely, at orbital position angles equal to 90° and 270°) 
This emphasizes the strength of polarimetry for extrasolar planet detection and 
characterization. Incidentally, is smallest for the inclination angles and orbital 
position angles where it is the most difficult or even impossible to directly observe 
the planet, i.e. at large inclination angles and orbital position angles equal to 0°, 
180°, or 360°, when the planet is close to, or even in front of or behind its star. 

4.2 Clear and cloudy planets with wavelength dependent 
surface albedos 

4.2.1 Wavelength dependence 

Figure [8] shows the wavelength dependence of the total flux, F, and degree of 
polarization, Ps, of starlight that is reflected by planets that are completely covered 
by, respectively, ocean and deciduous forest, with atmospheres that are either clear, 
i.e. cloudless, or cloudy, i.e. that contain a homogeneous cloud layer. The cloud 
and the cloud particles have been described in Sect. l3.1l For comparison, we have 
also plotted F and Ps of the clear white and black planets discussed in Sect. 14.1 1 
The planetary phase angle, a, is 90°. First, we will discuss F and Pj of the planets 
with the clear atmospheres, and then those of the cloudy planets. 

The surface albedo of the clear planet that is covered by ocean equals zero, 
regardless of wavelength (see Sect.[3]2]l- The differences between F and Ps of the 
clear black planet and F and Ps of the clear, ocean covered planet (see Fig. [Sj, 
are thus due to the specular reflecting interface between the model atmosphere and 
the ocean, which increases the total amount of light that is reflected back towards 
the observer. In Fig. [8^, the specular reflection increases F in the continuum with 
about 10% at the short wavelengths and with about 20% at the long wavelengths 
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(where more incoming starlight reaches the surface, because of the smaller atmo- 
spheric optical thickness). Although very difficult to see in Fig.[8ji, the influence of 
the specular reflection on F is small within the gaseous absorption bands, because 
at those wavelengths, little light reaches the surface and after a reflection there, the 
top of the atmosphere again. 

As can be seen in Fig.[8}3, the specular reflection decreases in the continuum 
with a few percentage points, because with the specular reflection, a fraction of 
the light that is incident on the surface is reflected back towards the atmosphere, 
adding mainly to the unpolarized flux (at least at this phase angle). In gaseous 
absorption bands, the influence of the specular reflection on is smaller than in 
the continuum, because at these wavelengths less light reaches the surface and after 
a reflection there, the top of the atmosphere again. 

The wavelength variation of the continuum F and P^ pertaining to the clear, 
forest-covered planet reflects the wavelength variation of the surface albedo (see 
Fig. O except at the shortest wavelengths. There, F and P^ are mainly deter- 
mined by the light that has been scattered by the gaseous molecules within the 
planet's atmosphere, because (1) the atmospheric scattering optical thickness in- 
creases with decreasing wavelength, and (2) the surface albedo is only about 0.05 
at the shortest wavelengths (see Fig. [3}. At longer wavelengths, the character- 
istic reflection by chlorophyll (around 0.54 yum) and the red edge (longwards of 
0.7 jum) can easily be recognized both in F and in P^. The red edge in flux spec- 
tra of the Earth has been detected by instruments onboard the Galileo mission on 
its way to Jupiter (Sagan et al. 1993), and from the ground it has been obs erved 
and in some cases m odeled by d ifi'erent research groups (s ee e.g. Iffamdani et al.l 
20061: iMontanes-Rodrfguez et al]|2006l ; iTinetti etani2006bl R ISeager et al.ll2005l : 



Woolf et al.ll2002l : [Arnold et al.ll2002h in spectra of Earthshine, the sunlight that 



has first been reflected by the Earth and then by the moon, and that can be ob- 
served on the moon's nightside. Interestingly, the reflection by chlorophyll leaves 
a much stronger signature in P^ than in F, because in this wavelength region P^ 
appears to be very sensitive to small changes in As, as can also be seen in Fig.|4j3. 

Adding a cloud layer to the atmosphere of a planet covered with either veg- 
etation or ocean, increases F across the whole wavelength interval (see Fig. [8^). 
A discussion on the effects of d ifferent types of clou ds on flux spectra of light re- 
flected by exoplanets is given bv lTinetti et all j2006bl lal). Our simulations show that 
although the cloud layers of the two cloudy planets have a large optical thickness 
(i.e. 10 at /I = 0.55 yum, as described in Sect. 13. lb both cloudy planets in Fig. [8^ 
are darker than the white planet with the clear atmosphere (the flux of which is 
also plotted in Fig. [8^). The cloud particles themselves are only slightly absorb- 
ing (see Sect. 13. U . Apparently, on the cloudy planets, a significant amount of 
incoming starlight is diffusely transmitted through the cloud layer (through mul- 
tiple scattering of light) and then absorbed by the planetary surface (the albedos 
of the ocean and the forest are smaller than 1.0). Thus, even with an optically 
thick cloud, the albedo of the planetary surface still influences the light that is 
reflected by the planets, and approximating clouds by isotropically or anisotrop- 
ically reflecting surfaces, without regard for what is underneath, as is sometimes 
done (see e.g. Montanes-Rodrfguez et al. 2006; Woolf et al. 2002) is not appropri- 
ate. Assu ming a dark surface b eneath scattering clouds with non-negligible optical 
thickness jTinetti et al.ll2006bl laD will lead to too dark planets. The influence of the 
surface albedo is in particularly clear for the cloudy planet that is covered with 
vegetation: longwards of 0.7 yum, the continuum flux of this planet still shows 
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the vegetation's red edge. The visibihty of the red edge through optically thick 
clouds strengthens the detecta bility of surface bio signatures in the visible wave- 
length range as discussed by jTinetti et alj|2006"bh . whose numerical simulations 
showed that, averaged over the daily time scale, Earth's land vegetation would be 
visible in disk-averaged spectra, even with cloud cover, and even without account- 
ing for the red edge below the optically thick clouds. Note that the vegetation's 
albedo signature due to chlorofyll, around 0.54 fjm, also shows up in Fig. [8^, but 
hardly distinguishable. 

The degree of polarization, f s, of the cloudy planets is low compared to that of 
planets with clear atmospheres, except at short wavelengths. The reasons for the 
low degree of polarization of the cloudy planets are (1) the cloud particles strongly 
increase the amount of multiple scattering of light within the atmosphere, which 
decreases the degree of polarization, (2) the degree of polarization of light that is 
singly scattered by the cloud particles is generally lower than that of light singly 
scattered by gaseous molecules, especially at single scattering angles around 90° 
(see Fig.|2j5), and (3) the direction of polarization of light singly scattered by the 
cloud particles is opposite to that of light singly scattered by gaseous molecules 
(see Fig.[2]5). Thanks to the latter fact, the continuum of the cloudy planets is 
negative (i.e. the direction of polarization is perpendicular to the terminator) at the 
longest wavelengths (about -0.03 or 3 % for /I > 0.73 yum in Fig. [8J5). At these 
wavelengths, the atmospheric molecular scattering optical thickness is negligible 
compared to the optical thickness of the cloud layer, and therefore almost all of the 
reflected light has been scattered by cloud particles. 

Unlike in the flux spectra, the albedo of the surface below a cloudy atmosphere 
leaves almost no trace in of the reflected light: the red edge of the vegetation 
hardly leads to a difference between P^ of the cloudy planets (Fig.[8j)). In par- 
ticular, at 1.0 jim, Ps of the cloudy, vegetation-covered planet is -0.030 (-3.0%), 
while f s of the cloudy, ocean-covered planet is -0.026 (-2.6%). The reason of 
the insensitivity of P^ of these two cloudy planets to the surface albedo is that the 
light reflected by the surfaces in our models adds mainly unpolarized light to the 
atmosphere, in a wavelength region where P^ is already very low because of the 
clouds. 

The cloud layer has interesting effects on the strengths of the absorption bands 
of O2 and H2O both in F and in P^: because the cloud particles scatter light very 
efficiently, their presence strongly influences the average pathlength of a photon 
within the planetary atmosphere. At wavelengths where light is absorbed by atmo- 
spheric gases, clouds thus strongly change the fraction of light that is absorbed, and 
with that the strength of the absorption band. These are well-known effects in Earth 
remote-sensing; in particular the O2 A-band is used to der ive e.g. cloud top alti- 
tudes and/or cl oud coverage within a ground pixel (see e.g. Kuze & Chance 1994: 
iFischer & Gra ssl 1991; Fischer et al. 1991; Saiedv et al. 1967; Stam et al. 2000^), 
because oxygen is well-mixed within the Earth's atmosphere. In general, clouds 
will decrease the relative depth (i.e. with respect to the continuum) of absorp- 
tion bands in reflected flux spectra (see Fig. [8^), because they shield the absorbing 
gases that are below them. However, because of the multiple scattering within the 
clouds, the absorption bands will be deeper than expected when using a reflecting 
surface to mimic the clouds. For example, the discrepancy between absorption 
band depths in Earth-shine flux o bservations and model simulations as shown by 
iMontafies-Rodriguez et alj ^20061) . with the observation yielding e.g. a deeper O2- 
A band than the model can fit, can be due to neglecting (multiple) scattering within 
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the clouds, as lMontanes-Rodriguez et ^ ( 12006!) themselves also remark. 

Another source for diiferences between absorption band depths in observed 
and modeled flux spectra could be that when modeling albedo and/or flux spec- 
tra, th e state of polarization of the light is usually neglected. IStam & HovenieJ 
l l2005l) showed for Jupiter-like extrasolar planets that neglecting polarization can 
lead to errors of up to 10 % in calculated geometric albedos, and that in particular 
the depths of absorption bands are affected, because the error in the continuum is 
usually much larger than the error in the deepest part of the absorption band. In 
Sect.|5] we will show that neglecting polarization does not significantly change the 
depth of the absorption bands in the flux spectra of these Earth-like model planets. 

In the polarization spectra (Fig. [HJ)), the effects of clouds on the strength of 
the absorption band features are more complicated than in flux spectra, because 
the absorption not only changes the amount of multiple scattering that takes place 
in the atmosphere, but also the altitude where the reflected light has obtained its 
state of polarization; in an inhomogeneous atmosphere, like a cloudy one, different 
particles at different altitudes can leave the light they scatter in differe nt states of 
polar ization (this also affects the scattered fluxes, but less so) (see e.g. IStam et al] 
Il999h . In Fig. [8J5, P, slightly increases within absorption bands in wavelength 
regions where the continuum is positive (A < 0.73 yum), whereas Ps slightly 
decreases (in absolute sense) within absorption bands in wavelength regions where 
the continuum P^ is negative (A > 0.73 ^m). In these cloudy model atmospheres, 
both the increase and the decrease (in absolute sense) of Ps in absorption bands 
are due to (1) a decrease of multiple scattering (which also takes place in purely 
gaseous model atmospheres), and (2) an increase of the relative amount of photons 
that are scattered by gaseous molecules instead of by cloud particles, since the 
latter are located in the lower atmospheric layers. 

The change of the strength of an absorption band in F or P^ due to the presence 
of a cloud layer, depends strongly on the altitude of the cloud layer, its optical 
thickness, the cloud coverage , the mixing ratio and the vertical distribution of 
the absorbing gas. For example, in both the flux and the polarization spectra of 
the cloudy planets in Fig. [8] the absorption bands of H2O are weak compared to 
the same bands in the spectra of the cloudless planets, because most of the H2O 
is located below the clouds. The absorption bands of O2 are also weaker for the 
cloudy planets than for the cloudless planets, although the influence of the clouds 
on these absorption bands is less strong than on the bands of H2O, simply because 
O2 is well-mixed throughout the atmosphere, and thus not primarily located below 
the clouds, like H2O. 

When a terrestrial type extrasolar planet will be discovered, it will of course 
be extremely interesting to try to identify oxygen in the planet's atmosphere, and 
in particular to determine the oxygen mixing ratio from absorption bands such as 
the O2 A-band. As discussed above, the depth of such an absorption band will 
depend not only on the absorber's mixing ratio, but also on the cloud cover. As an 
example. Fig. |9] shows the influence of the cloud top altitude and the O2 mixing 
ratio on the depth of the O2-A absorption band, both for the flux and the degree of 
polarization of the reflected starlight, for a = 90°. The figure shows the results of 
numerical calculations for model planets covered by ocean and with cloud layers 
of optical thickness 10 (at 0.55 /jm) placed with their tops at, respectively, 802 hPa, 
628 hPa, and 487 hPa, and with oxygen mixing ratios of, respectively, 11 %, 21 %, 
and 31 %. From Fig.|9ji, it is clear that the continuum flux is independent of the 
O2 mixing ratio, as it should be, and that it is virtually independent of the cloud 
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top altitude (note the vertical scale). The latter is easily understood by realizing 
that at these wavelengths, the gaseous molecular scattering optical thickness above 
the different cloud layers is negligible compared to the scattering optical thickness 
of the cloud layers themselves. Not surprisingly, the flux in the absorption band 
increases significantly with a decreasing O2 mixing ratio and with increasing cloud 
top altitude (i.e. decreasing cloud top pressure). Concluding, the depth of the O2 
A~band in planetary flux spectra depends on the O2 mixing ratio as well as on the 
cloud top altitude, and it can thus not be used for deriving the O2 mixing ratio 
if the cloud top altitude is unknown or vice versa. Figure |9j5 shows the relation 
between in the continuum and in the absorption band, and the cloud top altitude 
and oxygen mixing ratio. Apparently, the continuum depends significantly on 
the cloud top altitude, because the degree of polarization is very sensitive to even 
small amounts of gaseous molecules, and depends only very slightly on the oxygen 
mixing ratio. The degree of polarization in the absorption band depends on the 
oxygen mixing ratio as well as on the cloud top altitude. Concluding, the strength 
of the O2 A band in planetary polarization spectra can be used to derive both the 
cloud top altitude and the oxygen mixing ratio. The influence of e.g. broken cloud 
layers, and clouds at difl"erent altitudes on F and P^, will be subject for later studies. 

4.2.2 Phase angle dependence 

Figures [To] and [TT] show the phase angle dependence of F and P^ of the starlight 
that is reflected by the clear and cloudy, ocean and vegetation-covered planets 
appearing in the previous section (see Sect. 14.2.11 and Fig. [8]l at /i = 0.44 fim 
(Fig. llOl l and at A = 0.87 yum (Fig. lilt , respectively. Remember that at phase 
angle a = 0°, the fluxes plotted in Figs. 1 10b andlllb are just the planets' geometric 
albedos, Aq, at those wavelengths. For the clear, ocean-covered planet, Aq is 0.15 
at A = 0.44 yum, and 0.014 at 0.87 yum. For the clear, vegetation-covered planet, 
Aq is 0.16 at 0.44 fim, and 0.37 at 0.87 yum. For the cloudy, ocean-covered planet, 
Aq is 0.49 at 0.44 //m, and 0.52 at 0.87 ^m, and for the cloudy, vegetation-covered 
planet, Aq is 0.49 at 0.44 fim (like for the cloudy, ocean-covered planet), and 0.60 
at 0.87 yUm. 

For the two cloudy planets, the reflected flux is not a smoothly decreasing 
function of phase angle, but instead shows some bumpy features near a = 5° 
and 30° for A = 0.44 yum (Fig. [Toji), and near 10° and 35° for A = 0.87 yum 
(Fig. II lb). These angular features trace back to features in the single scattering 
phase function of the cloud particles (at scattering angles, 0, larger than 140°) (see 
Fig. |2ji). Because at A = 0.44 yum, the gaseous molecules above the cloud layer 
scatter light more efliciently than at A = 0.87 yum, the angular features in F are 
more subdued at 0.44 yum than at 0.87 yum. Interestingly, at phase angles near 60° 
and for A = 0.87 yum (Fig.ll lb), the cloudy, ocean-covered planet is about as bright 
as the clear, vegetation-covered planet. Apart from the high surface albedo of the 
vegetation at this wavelength (Fig. [3]l, this can partly be attributed to the single 
scattering phase function of the cloud particles, too, because this function has a 
broad minimum between = 90° and = 120° (see Fig.|2^). Note further that 
only at the largest phase angles, i.e. a > 120°, the cloudy planets are brighter than 
the clear planet with the surface albedo equal to 1.0. This is in particularly obvious 
at A = 0.87 yum (Fig. I lib), and is due to the strong forward scattering peak in the 
single scattering phase functions of the cloud particles (Fig. |2^). At these large 
phase angles, the light that is reflected by the cloudy planets has not reached the 
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surface, and the reflected flux is therefore independent of the surface albedo (see 
Fig.EJi). 

The phase angle dependence of the degree of linear polarization, P^, of the light 
reflected by the two cloudy planets (Figs.llOb andll lb) shows, for a < 40°, angular 
features that are due to angular features in the single scattering polarization phase 
function of the cloud particles (see Fig.[2j3), just like the phase angle dependence 
of the reflected flux (Figs.[TOji and[nii)- At A = 0.44 /jm, Ps is determined not only 
by light scattered by the cloud particles, but also by light scattered by the gaseous 
molecules, whereas at A = 0.87 yum, P^ is predominantly determined by the cloud 
particles. As a result, al A = 0.87 fim, the characteristic polarization signatures of 
light scattered by the cloud particles (i.e. the angular features at ff < 40° and also 
the negative values for 60° < a < 160°) are much stronger than at A = 0.44 yum 
(taking into account the wavelength dependence of the single scattering features 
themselves, cf. Fig.|2j. 

The angular features around a = 32° for A = 0.44 fim, and around a = 25° for 
A = 0.87 fim, pertain to the so-called primary rainbow feature. AtA = 0.44 /im, the 
degree of polarization of this feature is about 0. 1 (10%) and alA = 0.87 fim, about 
0.06 (6% ). This is muc h lower than the 20% predicted, for a completely cloudy 
Earth, bv lBailevI ( l2007h . This discrepancy is probably due to the fact that lSailevI 
io^ arrives at his disk-integ rated value by integratin g local observations by the 
POLDER-satellite instrument jPescharnps etal.lil994l) . while our multiple scat- 
tering calculations and disk-integration method take into account the variation of 
illumination and viewing angles across the planetary disk. In addition, dilferences 
in single-scattering properties of cloud particles (which depend on the particle size 
distribution and composition) and the optical thickness of the clouds will influence 
the strength of the feature, although probably not more than a few percent. 

At neither of the two wavelengths, P^ of the cloudy planets is sensitive to the 
albedo of the surface below the clouds (Figs.fTOb andll lb): only at 0.87 /im and 
for a < 40°, Pg of the cloudy, ocean-covered planet is at most 0.01 (1%) larger 
than Ps of the cloudy, vegetation-covered planet, mainly because of the darkness 
of the ocean. 

For the clear planets with the Lambertian reflecting surfaces, the phase angle 
dependence of Ps at both wavelengths (Figs.fTOb andll lb) is similar to that shown 
in Fig. [5J3, taking into account the differences in surface albedo. Compared with 
the clear, black planet, the specular reflecting surface of the clear, ocean-covered 
planet decreases Ps at all phase angles (with at most 0.04 at o- = 90° and A = 
0.87 fim), except at phase angles larger than about 130° at 0.44 ftm, and about 
145° at 0.87 //m. 

4.3 Clear and cloudy planets with horizontal inhomo- 
geneities 

4.3.1 Wavelength dependence 

In Sects. |4n and |42l we have presented numerically simulated fluxes and de- 
grees of polarization of light reflected by horizontally homogeneous planets. In 
this section, we will show results for quasi horizontally inhomogeneous planets, 
using Eq.|15|and the flux vectors calculated for the clear and cloudy, ocean- and 
vegetation-covered planets presented in Sect l4.2l and Figs.[8l[T0l and[TT] The re- 
sulting spectra can be thought of as the Earth observed as if it were an exoplanet, 
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and using an integration time of at least a day. Since we use a weighted sum of 
liomogeneous planets, our spectra might differ from those obtained using a model 
planet covered by continents and oceans, even if the coverage fractions are the 
same. Our spectra will, however, give a good estimate of what might be expected. 

Although of course endless combinations of these flux vectors can be made, we 
will limit ourselves here to Earth-like combinations, and leave other combinations 
and retrieval algorithms for subsequent studies. In addition, including more types 
of surface coverage, different types of clouds, and e.g. different cloud coverages for 
different surface types, would add d etails to the m odelling results that are beyond 
the scope of this article. See e.g. Ixinetti et for examples of flux 

spectra for different cloud types and surface coverages. 

Figure [12] shows the flux and degree of polarization of light reflected by an ex- 
oplanet that has, like the Earth, 70 % of its surface covered by a specular reflecting 
ocean and 30 % by deciduous forest. The cloud coverage ranges from 0.0 (a clear 
atmosphere) to 1.0 (a completely overcast sky) in steps of 0.2. Recall that the mean 
global cloud coverage of the Earth is about 0.67 jRossow et al.l 19961) . Note that to 
simulate F and of the cloudy fractions of the planets, we combined the cloudy, 
ocean-covered planets with the cloudy, forest-covered planets. In Eq.[T5] we thus 
used A' = 4. 

The main, not surprising difference between the flux of the clear, quasi hori- 
zontally inhomogeneous planet from Fig. 112b . and that of the clear, horizontally 
homogeneous forest-covered planet in Fig.[8ji, is that the red edge (for A > 0.7 yum) 
that is characteristic for reflection by vegetation, is much less strong when 70 % of 
the planet is covered by ocean than when the whole planet is covered by vegeta- 
tion (the red edge continuum is approximately 0.035 for the inhomogeneous planet 
and 0.11 for the homogeneous one). The (black) ocean all but removes the local 
maximum in F at green wavelengths (between 0.5 and 0.6 yum) which is due to 
chlorophyll in the vegetation. In (Fig. 112b). the ocean somewhat changes the 
spectral shape of the red edge feature, and decreases its depth by about 0.04 (at 
A = 0.87 yum) when compared to Fig.[8l5. The ocean signficantly changes the spec- 
tral feature in that is due to chlorofyll: for the clear, completely forest-covered 
planet (Fig.[8]5), the minimum across this feature is 0.38 while it is 0.63 for the 
clear planet in Fig. 1 12b. Finally, in P,, the gaseous absorption bands are stronger 
for the clear planet covered by ocean and forest, than for the clear, forest-covered 
planet. 

From Fig. 112b. it is clear that F in the continuum is very sensitive to the cloud 
coverage. The continuum is also sensitive to the cloud coverage, but this sensi- 
tivity decreases with increasing cloud coverage, in particularly at the longer wave- 
lengths. 

To get more insight into the sensitivity of the shape of the flux and polarization 
spectra to the surface and the cloud coverage, we have plotted in Fig. [13] F and 
Ps in the near-infrared continuum (A = 0.87 yum) against F and Ps in the blue 
continuum (A = 0.35 fim) for planets with surface coverage ratios ranging from 
0.0 (100 % forest) to 1.0 (100 % ocean), in steps of 0.2. The cloud coverage 
ranges from 0.0 (a clear planet) to 1.0 (a completely cloudy planet). Looking 
at the reflected fluxes (Fig. 113b). it can be seen that for a given cloud coverage, 
F in the blue is virtually independent of the surface coverage (ocean or forest), 
while in the near-infrared, the sensitivity of F to the surface coverage depends 
strongly on the cloud coverage: it is relatively large when the cloud coverage is 
small, and relatively small when the cloud coverage is large, as could be expected. 
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Fig. [T3h shows that our completely cloudy planets are somewhat brighter in the 
near-infrared than in the blue (assuming a wavelength independent stellar flux, or, 
after correcting observations for the incoming stellar flux). Looking at the degree 
of polarization of the reflected fluxes (Fig. 1 13b). it is clear that the larger the cloud 
coverage, the smaller the dependence of at both 0.35 yum and 0.87 yum on the 
surface coverage. For a given cloud coverage smaller than about 0.5, Ps can be seen 
to depend on the surface coverage. In particular, the larger the fraction of ocean 
on the planet, the larger in the near-infrared. In the blue, f s also increases with 
increasing fraction of ocean coverage, but only slightly so. 

It is important to remember that the precise location of the data points in 
Fig. [T3] and thus the retrieval opportunities, will depend on the physical param- 
eters of the model cloud, such as the cloud optical thickness, the microphysical 
properties of the cloud particles, and the altitude of the cloud. We will explore 
such dependencies in later studies. 

4.3.2 Phase angle dependence and diurnal variation 

As a horizontally inhomogeneous planet like the Earth rotates around its axis, the 
surface fraction of, for example, ocean that is turned towards a distant observer 
will vary during the day (except when the observer is located precisely above one 
of the planet's geographic poles). Assuming a planet with a surface that is covered 
only by land and water, the diurnal variation of the distribution of land and water 
across the part of the planetary disk that is illuminated and turned towards a dis- 
tant observer depends on many factors, such as the actual distribution of land and 
water across the planet, the sub-observer latitude, the planet's phase angle, and the 
location of the terminator (the division between day and night on the planet, that 
will depend on the obliquity of the planet and the time of year). The variation that 
can actually be observed, would, of course, also depend on the cloud cover. Cal- 
culated d iurnal vari ations of fluxes of Earth-like exoplanets have been presented 
bv ,Fordetal] ( l200lh . using a single scattering Monte C arlo radiative transfe r code 
and horizontally inhomogeneous model planets, and bv lTinetti et al.l j2006bh . 

Here, we will show a few examples of the diurnal variation of the flux and in 
particular the degree of polarization of a quasi horizontally inhomogeneous exo- 
planet with a longitudinal distribution of land and water similar to that along the 
Earth's equator (see Fig. I14t . For simplicity, we assume the distribution is lat- 
itude independent (we thus have a planet with "vertical stripes"). Furthermore, 
the distant observer is located in the planet's equatorial plane, and we assume 
the plan et's obliquity equ als zero. Our model plan ets are somewhat simpler than 
those of lFord et alj j200lh and lTinetti et al. I( l2006bl) , who use realistic ally horizon- 
tally inh omogeneous planets, but we do use multiple scattering, like lTinetti et al.l 
( l2006bl) . and polarization. 

In Fig. [TS] we have plotted F and of a cloudless planet seen under phase 
angles of 50° (when more than half of the illuminated part of the planetary disk 
is visible), 90° (quadrature, when half of the illuminated planetary disk is visible) 
and 130° (when less than half of the illuminated planetary disk is visible). We 
show curves for A = 0.44 yum, and for A = 0.87 nm. Figur e[T6]is similar to Fig. 1151 
except that here, like on an average Earth teossow et al.iri996h . two-thirds of the 
planet is covered by clouds. 

For the cloudless planet (Fig. |15l >, F and P^ show for all three phase angles 
much more variation in the near-infrared (0.87 yum) than in the blue (0.44 /jm). 
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This is not surprising, since in the blue, the vegetation and the ocean are both 
dark, and F and are thus mostly determined by the gaseous atmosphere. In 
the near-infrared, the contrast between the (red edge of the) vegetation and the 
ocean shows up clearly in both F and P^. A similar wavelength dependence of the 
variation of the reflected flux (at a = 90°) was found by Ford et al. (2001) and by 
iTinetti et alj i2p06b) (the latter have 0.50 to 0.55 /jm as the shortest wavelength 
region ). In the blue (0.44 /jm in our plots, and 0.45 fim in the plots bv lpord et al.l 
OOOlh). our reflecte d flux (for a = 90°) shows far less diurnal variation than that 
of Ford et al. I {20013. This is most likely due to the multiple scattering that we 
include in our calculations. In the infrared (0.87 /im in our plots, and 0.75 fim in 
the plots bylpp rd et al. (2001)), our maximum fluxes (for a = 90°) are larger than 
those of iFord et aU ([2001), namely about 0.06 versus about 0.03, which cannot 
be attributed to multiple scattering, because the atmospheric molecular scattering 
optical thickness at these wavelengths is very small. This difference in maximum 
flux is most probably due to the difference in surface coverage (vegetation versus 
Sahara sand and ocean) and hence in surface albedo at these wavelengths. 

The diurnal variation in the degree of polarization, P^, has not been studied 
before. In the blue, P^ is strongly determined by Rayleigh scattering, and thus 
depends strongly on the planetary phase angle. Here, the average P^ for a = 50° 
is very similar to that for a = 130°, because of the symmetry of the degree of 
polarization of light that is singly scattered by gaseous molecules with the single 
scattering angle (see Fig. [2J3 and Fig. llOb). The variation in Ps is largest (in ab- 
solute sense) for a = 90°, and can be traced back to the difference, at o- = 90°, 
between the curves for the vegetation-covered and the ocean-covered planets in 
Fig. 1 10b . Furthermore according to the difference between these two curves in 
Fig.lIOb. the diurnal variation for a = 50° should be significantly larger than that 
for a = 130°. However, as can be seen in Fig. [15] the diurnal variation for these 
two phase angles is similar, and even slightly larger for a = 130° than for a = 50°. 
The reason for this is that when a = 130°, the illuminated and observable part of 
the planet is smaller than for a = 50°, and will therefore have a more uniform 
surface coverage (either ocean or vegetation) at any given time. 

At red wavelengths (lower panels in Fig. |15t , where the largest contribution 
to the reflected light is from the surface, the diurnal variation of the degree of 
polarization, P^, depends strongly on the planetary phase angle; the variation in Ps 
clearly increases with the phase angle. 

Next, we discuss the influence of clouds on the diurnal variation of F and 
Ps- As can be seen in Fig.[T6l and as can be expected from Fig. (TO] in the blue 
(0.44 yCim), our clouds, which cover two-thirds of the planet, smother virtually any 
diurnal variation in F and in P^. In the near-infrared (0.87 fim), the flux F retains 
some of the diurnal variation, because at this long wavelength the vegetation has 
a much higher albedo than the ocean (cf. Fig. II It. The fluxes we c alculate for 
the cl oudy planets (and a = 90°) are very similar to those calculated bv lFord et al.l 
j2001l) , because to describ e the light reflected by the cloudy regions on their model 
planets, iFord et al.U200ll) apparently do not rely on the single scattering approach 
they use for the cloud-free planets. Instead they adopt the reflected fluxes presented 
by [Hovenier & Hage ( 1989) for cloudy planet models, the calculation of which 
does fully include multiple scattering (though no vegetated surfaces underneath 
the clouds). Qualitatively, our variation in F is very similar to that shown by 
iTinetti et alj l2006b), although the normalization is different. 

The degree of polarization in the near-infrared shows no measurable variation; 
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the higher albedo of the vegetation only adds more unpolarized light to the total 
reflected light. All clouds on our model planet are optically thick (b'^ = 10 at 
0.55 ^m, see Sect.[3}. Optically thinner clouds, that depolarize the light less, would 
probably leave more diurnal variation in P^, as would e.g. polarizing vegetation 
surfaces (Wolstencroft et al. 2007). 

With our method for calculating the light reflected by horizontally inhomoge- 
neous planets (i.e. using Eg. 115b. light reflected by e.g. the cloud layer effectively 
comes from all locations on the illuminated and observed part of the planet, in- 
stead of from a limited number of locations spread over the planet, which would 
be a more realistic cloud deck. We might therefore overestimate the contribution of 
the cloud layer to the total reflected signal. The difference between the signals of 
our quasi inhomogeneous planets and those of more realistically inhomogeneous 
planets require more investigation, also with the view on the retrieval of planetary 
characteristics. Such an investigation should also address the region-to-region vari- 
ations in time averaged cloud coverage and cloud optical thickness, that are found 
on Earth when cloud variatio ns on time scales less than a month (i.e. weather) are 
removed (for cloud data, see lRossow & Dueiiasll200 4^. Although the actual spa- 
tial distribution of cloud properties across an extrasolar planet might not resemble 
that of the Earth, relations between land and cloud coverage. Regarding flux cal- 
culations, examp les of region-to-region variations in cloud coverage are given by 
lFordetalJ ( l200ll) . 

5 Summary and conclusion 

We have presented numerical simulations of the flux and state (degree and direc- 
tion) of polarization of starlight reflected by various types of Earth-like extrasolar 
planets as functions of the wavelength and as functions of the planetary phase 
angle. S pectral fluxes of Earth-like extrasolar planets have been presented before 
(see, e.g. lTinetti et alj2006a c; Montanes-Rodriguez et al.l2006l . [Tumbull et alj2006l 
and references therein), the spectral variation of the degree of polarization of such 
planets not. 

Our results clearly show that the light reflected by an Earth-like exoplanet can 
be highly polarized, with the degree and direction of polarization depending on the 
physical characteristics of the planetary atmosphere and surface, on the illumina- 
tion and viewing geometries, and on the wavelength. Polarimetry thus appears to 
be a useful tool to distinguish the (polarized) stellar light that has been reflected 
by an Earth-like extrasolar planet from the (unpolarized) direct stellar light, and 
hence to detect Earth-like extrasolar planets. 

Our results also show that for given physical properties of the Earth-like ex- 
trasolar planet, the degree of polarization of the reflected light has similar spectral 
features as the flux of this light: there are features related to the surface albedo 
of the planet and high spectral resolution features that are due to absorption of 
light by atmospheric gases. The latter features, although due to different gases, 
were also present in cal culated polarizati on spectra of starlight reflected by ex- 
trasolar gaseous planets l lStam et ^l2004h . For Earth-like extrasolar planets, the 
occurence of spectral features due to gaseous absorption in polarization spectra is 
especially interesting because these features are conserved upon transmission of 
the light through the Earth's atmosphere even if it has a similar composition as 
the extrasolar planet's atmosphere (although the number of available photons in 
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the absorption bands will be greatly reduced upon traveling through the Earth's 
atmosphere if it contains the same absorbing gases as the extrasolar planet's atmo- 
sphere). 

Interestingly, the degree of polarization of reflected light appears to have a dif- 
ferent sensitivity to the structure and composition of the planetary atmosphere and 
the albedo of the underlying surface than the flux of the reflected light has. Com- 
pared to using only spectroscopy, polarimetry could thus provide additional and 
different information on the structure and composition of an extrasolar planet's at- 
mosphere and surface. In particular, we have shown that the degree of polarization 
at continuum wavelengths around the O2 A-band (around 0.76 nm) is sensitive to 
the cloud top altitude, whereas the continuum flux in this wavelength region is vir- 
tually insensitive to the cloud top altitude. The degree of polarization could thus 
be used for cloud top altitude determination, a method th at is also applied to e.g . 
observations of the Earth observing POLDER instrument jPeschamps et al J 19941) 
and its successors. We have also shown that the depth of the O2 A-band in a po- 
larization spectrum is sensitive to the oxygen mixing ratio, whereas the flux in the 
O2 A-band is sensitive to both the cloud top altitude and the oxygen mixing ratio. 
Polarimetry could thus help to disentangle crucial information on gaseous mixing 
ratios and clouds from the sparse, spatially unresolved data that will be available 
for exoplanets. 

Because of the strengths for detecting and characterizing Earth-like exoplan- 
ets, as mentioned above, polarimetry is a technique used in SEE-COAST (the Su- 
per Earths Explorer ~ Coronographic Off- Axis Telescope), a space-based telescope 
for the detection and the characteriz ation of gaseous exopl anets and large rocky ex- 
oplanets, so called 'Super-Earths', jSchneider et alj2006l) , that has been proposed 
to ESA in response to its 2007 Cosmic Vision call. 

Our numerical simulations should not only be useful for researchers that are 
interested in designing and building (spectro)polarimeters for the detection and 
characterization of Earth-like exoplanets, and for researchers interested in measur- 
ing and analyzing the state of polarization of extrasolar planets. Indeed, due to 
their optics, spectrometers tend to be sensitive to the state of polarization of the in- 
coming light (unless carefully corrected for): the measured fluxes thus depend on 
the deg ree and direction of polarization of the observed light (see e.g. IStam et al] 
l2000bl) . In case the polarization sensitivity of a spectrometer is known, e.g. be- 
cause it has been measured during calibration, our numerical simulations can help 
to estimate the error in the measured flux. To actually correct measured fluxes 
for the polarization sensitivity of an instrument, one has to know the instrument's 
polarization sensitivity and one has to measure the state of polarization of the in- 
coming light, because as our simulations show, this state of polarization varies 
with the illumination and viewing geometries and with the physical characteristics 
of the planet. Because the state of polarization is wavelength dependent and shows 
high spectral resolution features similar to those in the total flux, the state of po- 
larization should be measured with the same spectral resolution as the total flux 
if one is interested in accurately measuring the depth and shape of e.g. gaseous 
absorption bands. 

Finally, IStam & HovenieJ j2005t) showed that neglecting polarization when 
calculating disk-integrated total flux spectra of gaseous planets induces errors of 
several percent across the continuum and gaseous absorption bands. The reason for 
these errors is that the flux of scattered light depends on the state of po larization of 
the incident light (see also lMishchenko et~aill 19941 : iLacis et alj|l998l) . and while 
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the incident stellar light can be assumed to be unpolarized jKempetalJ 19871) , It IS 
usually polarized upon the first scattering within the atmosphere. Hence, as long as 
only single scattering processes are taken into account, polarization does not influ- 
ence the scattered total flux, and can thus safely be ignored if one is only interested 
in total fluxes. When multiple scattering has to be accounted for, however, neglect- 
ing polarization and treating light as a scalar instead of as a vector, does influence 
the scattered total flux. The fluxes presented in this paper have all been calculated 
taking polarization fully into account. To investigate the influence of neglecting 
polarization on the fluxes of Earth-like extrasolar planets, we performed a number 
of calculations without polarization. 

Figure [T7] shows total fluxes calculated with and without polarization, for a 
cloudfree and a cloudy planet, as functions of the wavelength at quadrature, and at 
A = 0.35 fjm as functions of the phase angle. For a cloudfree atmosphere with the 
surface albedo of 0.4, the maximum error is about 4%, at a wavelength of 0.35 /im. 
For the cloudy atmosphere with the vegetation underneath, the maximum error is 
about 2.5 %, also near 0.35 ^m. At other phase angles (Fig.ll7b), the relative errors 
have similar sizes, although they vanish near a = 50° and 120°. The errors have a 
similar ph ase angle dependence, bu t are smaller than those presented for gaseous 
planets bv lStam & Hovenie3 ( l2005h . The errors are smaller because our Earth-like 
atmosphere is optically much thinner than that of a gaseous planet, and hence less 
multiple scattering takes place. That for an Earth-like planet (with a total molecu- 
lar scattering optical thickness of about 0.1 at /t = 0.5 5 //m, see TablefTJ, the errors 
will be smaller is also clear from Fig. 5 o f IStam & H ovenier (2005), which shows 
the errors due to neglecting polarization in planetary albedos as functions of the 
atmospheric molecular scattering optical thickness. The decrease of the error with 
decreasing molecular scattering optical thickness also explains the decrease of the 
error with wavelength in Fig. [17] and the smaller error for the cloudy atmosphere. 
Concluding, for Earth-like model atmospheres, neglecting polarization when cal- 
culating fluxes leads to maximum errors of a few percent in the continuum. 
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star 




towards the observer 



Figure 1 : Top-view sketch of the geometries involved in observing extrasolar planets: 
D is the distance between the star and the planet, R is the radius of the star, d is the 
distance between the planetary system and the observer, r is the radius of the planet, 
and a is the planetary phase angle. We assume that D » R and that d » r. 
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Figure 2: The phase function (scattering mattix element Fn) (on the left) and the de- 
gree of hnear polarization (-F21 /F\i) (on the right) as functions of the scattering angle 
©, for light singly scattered by gaseous molecules and cloud droplets. The wavelength 
of the hght is 0.44 jum (dashed hne, only for the cloud droplets), 0.55 jum (solid lines), 
and 0.87 jum (dotted line, only for the cloud droplets). 
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Figure 3: The measured albedos of three types of common terrestrial vegetation as 
functions of the wavelength: conifers (dashed Une), deciduous forest (solid line), and 
grass (dotted line) (data from the ASTER spectral library). For our model planets, we 
use the albedo of deciduous forest to represent vegetated surfaces. 
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Figure 4: The flux F (left) and the degree of linear polarization Ps (right) of starlight 
reflected by model planets with clear atmospheres and isotropically reflecting, com- 
pletely depolarizing surfaces as functions of the wavelength, for various values of the 
(wavelength independent) surface albedo: 0.0, 0.1, 0.2, 0.4, 0.8, and 1.0. The planetary 
phase angle a is 90°. 
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Figure 5: The flux F (left) and flie degree of linear polarization Ps (right) of starlight 
reflected by the model planets with As - 0.0, 0.4, and 1 .0 of Fig. |4]as functions of the 
phase angle a. F and Ps have been plotted for two wavelengths: 0.44 fim (solid lines) 
and 0.87 /im (dashed lines). 
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Figure 6: The fluxes of Fig.|5t, normalized to 1.0 at a = 0°. For A = 0.87 fim (dashed 
lines), the curves for As = 1 .0 and As = 0.4 are virtually indistinguishable, and follow 
the theoretical (normalized) curve e xpected for a Lambertian reflecting sphere (see 
van de Hulst 1980t Stam et al. 20061) . i.e. F(a) - - (sin a + (n - a) cos a), with a in 



radians. 
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Figure 7: F (left) and (right) of stariight with A - 0.44 jim that is reflected by the 
model planet with = 0.4 as functions of the orbital position angle for the following 
orbital inclination angles i: 0° (dot-dashed lines), 30° (dashed Unes), 60° (dotted Unes), 
and 90° (soUd Unes). 
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Figure 8: The wavelength dependent flux, f , (left) and degree of polarization, P^, 
(right) of starlight that is reflected by clear and cloudy horizontally homogeneous 
model planets with surfaces covered by deciduous forest (thin solid lines), a specular 
reflecting ocean (thin dashed lines). Note that the lines pertaining to of the cloudy 
atmospheres are virtually indistinguishable from each other For comparison, we have 
also included the spectra of the clear model planets with surface albedos equal to 0.0 
and 1.0 (thick solid hnes), shown before in Fig.|4] The planetary phase angle, a, is 90°. 
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Figure 9; The flux F (left) and the degree of linear polarization (right) of starlight 
reflected by cloudy ocean planets, at /I = 0.762 yum, with absorption by O2 (along the 
x-axes) and without absorption by O2 (along the y-axes). The top of the cloud layer 
was located at 802 hPa (the 'low cloud'; the nominal altitude), 628 hPa (the 'middle 
cloud'), or 487 hPa (the 'high cloud'). The O2 mixing ratio was 1 1%, 21% (the nominal 
value), or 3 1%. The planetary phase angle a is 90°. 
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Figure 10: The flux F (left) and the degree of linear polarization (right) of starlight 
reflected by the model planets of Fig. [8] as functions of the phase angle a, for A - 
0.44 fjxn. The thin solid line pertains to the vegetation-covered planet, and the thin 
dashed line to the ocean-covered planet. The lines pertaining to the cloudy planets are 
virtually indistinguishable from each other. 
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Figure 11: Same as Fig. [10] except for A = 0.87 fim. The lines for pertaining to the 
cloudy planets are virtually indistinguishable from each other, except at phase angles 
between 20° and 40°. 
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Figure 12: Similar to Fig.|8] except for quasi horizontally inhomogeneous planets with 
70% of their surfaces covered by specular reflecting ocean and 30% by deciduous for- 
est. The cloud coverage ranges from 0.0 to 1.0, in steps of 0.2. The planetary phase 
angle, a, is 90°. 
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Figure 13: The relations between F (left) and (right) in the blue (A - 0.35 /zm) 
and the near-infrared {A = 0.87 fim) for planets with ratios between ocean and forest- 
surface coverage ranging from 0.0 (100 % covered by forest, indicated by the stars), 
to 0.0 (100 % covered by ocean, indicated by the black circles), with steps of 0.2. In 
addition, the cloud coverage ranges from 0.0 to 1 .0, in steps of 0. 1 . The planetary phase 
angle, a, is 90°. Note that we have included the line F(0.87 /im)=F(0.35 /im) in the 
graph on the left, and that at A - 0.87 fim, f ^ of the clear, 100 % ocean planet falls of 
the vertical figure scale; its value is 0.9. 
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Figure 14: The percentage of surface covered by water along the equator of the Earth. 
This curve is derived from data of the Clouds and the Earth's Radiant Energy System 
(CERES) instrument on NASA's Earth Observing System (EOS). The longitude of 0° 
corresponds to the meridian of Greenwich. 
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Figure 15: The flux F (left) and degree of polarization Ps (right) of starlight reflected by 
a cloudless quasi horizontally inhomogeneous planet as functions of the sub-observer 
longitude. The longitudinal distribution of land (covered by vegetation) and water 
(ocean) is as given in Fig. [14] The planetary phase angles are 50° (soUd lines), 90° 
(dashed lines) and 130° (dash-dotted lines). The upper two panels are for A = 0.44 fim, 
and the lower two panels for A = 0.87 /im. 
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Figure 16: Similar to Fig. [15] except for planets with a cloud coverage of 0.67 



dRossow et al. I ll996l) . The planetary phase angles are 50° (solid lines), 90° (dashed 



lines) and 130° (dash-dotted lines). The upper two panels are for A — 0.44 yum, and the 
lower two panels for A - 0.87 //m. 
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Figure 17: Total reflected fluxes as functions of the wavelength for a - 90° (left), 
and as functions of the phase angle for A - 0.35 /im (right). The solid lines have 
been calculated with polarization, and the dashed lines without polarization. The lower 
two lines are for a homogeneous, cloudfree planet with a surface albedo of 0.4, and 
the upper two lines are for a homogeneous, cloudy planet with a surface covered by 
vegetation. 
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Table 1 : The altitude z (in km), the pressure p (in hPa), the temperature T (in K), and 
the ozone (O3) and watervapor (H2O) m ixing ratios (in ppm, or p arts per million) at 
the 17 levels of the model atmospheres ( IMcClatchev et al.lll972h . Also given is the 
molecular scattering optical thickness b™^^ (at A - 0.55 ^m) of each of the 16 layers of 
the model atmospheres. The total molecular scattering optical thickness of the model 
atmosphere is 0.0975 (at A - 0.55 fim). Here, x(-y) stands for x ■ 10"' . 
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